1 INTRODUCTION 



The conceptual foundations of quantum mechanics continue to be a subject of 
debate and different interpretations of the theory still contend for acceptance. 
In this paper, we employ the the interpretation of quantum mechanics due to 
Lande[l,2] on spin theory and thereby to obtain a new and more generalized 
treatment of spin. 

In this work, we show how to derive the spin operators and eigenvectors from 
probability amplitudes. We obtain generalized forms of these quantities. To 
the best of our knowledge, both the method of derivation and the generalized 
expressions we present are new. The method shows that spin operators and 
vectors can be derived from an analysis which takes the probability amplitudes 
that characterize spin measurements as being primary, in the same way that 
the spherical harmonics and differential orbital angular momentum operators 
might be taken as being more fundamental than column vectors and matrices in 
describing orbital angular momentum. Thus the method we introduce does not 
require us to study the form of orbital angular momentum matrices and column 
vectors in order to construct spin vectors and operators by analogy. At the 
same time, this method places the treatment of spin on an even more analogous 
footing with that of orbital angular momentum than hitherto. 

This paper is constructed as follows. In Section 2, we give a brief review of 
the Lande approach to quantum mechanics, because it is likely to be unfamiliar 
to most readers. In Section 3 we use the Lande approach to investigate the 
transition from wave to matrix mechanics. In Section 4 we specialize the results 
of the previous sections to the case of spin- 1/2. We show in Subsection 4.1 how 
the matrix treatment of spin may be deduced from an approach that starts with 
probability amplitudes. In Subsection 4.3, we show how to obtain the Pauli spin 
matrices and their eigenvectors by this approach. We derive explicit formulas for 
the probability amplitudes characterizing spin-1/2 measurements in Subsection 
4.4. In Subsections 4.7, 4.8 and 4.9, we obtain the most general forms of the 
spin operators and give the eigenvectors of the generalized z component. In 
Subsection 4.10, we show that in the appropriate limit, the generalized operators 
and vectors reduce to the Pauli forms, and to those forms heretofore considered 
to be the most general . In Subsection 4.11, we discuss the connection between 
our results and the standard results. We conclude the paper after a Summary 
and Discussion in Section 5. 



2 REVIEW OF THE LANDE APPROACH 

2.1 Basic Results 

We begin by giving those results of the Lande approach which we shall need 
for our purpose. The contents of this section will be found set out at complete 
length in Refs.l and 2. 

According to Lande, we must start by accepting that indeterminism is inher- 
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ent in nature, and therefore in quantum measurement. Therefore upon repeated 
measurement, different results may follow from exactly the same preparation of 
a system. Statistical theory predicts the probabilities of obtaining the several 
different outcomes of any such measurement. 

We may briefly characterize quantum mechanics as dealing with material 
bodies composed of particles with coordinates and momenta q and p and with 
various observable quantities A,B,C,.. defined as functions of the q's and p's. 
Measurement forces the measured quantity to take one of the eigenvalues of the 
operator representing the measured quantity. The observable A can manifest 
many possible values A\, A 2 , A3, ... when it is measured. 

Suppose a system is in a state characterized by the eigenvalue Ai. A mea- 
surement of the quantity B is then made. Any one of the values of the eigenvalue 
spectrum of B can result. The probabilities of obtaining the different values of 
B are P(Af, Bj). Conversely, if we start with the system being in a state char- 
acterized by the eigenvalue Bj and make a measurement of A, then any of the 

values A\, A2, A3 can result. The probabilities of obtaining the different 

values of A are P{Bj\ Ai). The two sets of probabilities obey two-way symmetry. 
Hence 



P(A i ;B j ) = P{B j] A i ). 
The probabilities sum to unity: 



J2P(A i ;B j )=J2P(B j ;A i ) = l. 



(1) 



(2) 



To ensure reproducibility of results upon repetition of measurement, the prob- 
abilities satisfy 



P{Ai-A 6 



(3) 



We shall assume that both A and B have the same number of eigenvalues 
N; this can be shown to be always so [3]. We arrange the probabilities in a table 
to obtain the matrix 



[Pai 



( P{Ai;B{) P(A i; B 2 ) 
P(A 2] B 1 ) P(A 2 ;B 2 ) 

\ P(A N ;B X ) P{A N ;B 2 ) 



P(Ai;B N ) \ 
P(A 2 ;B N ) 

..P{A N -B N ) ) 



(4) 



Each row sums to unity, and, due to the symmetry of the probabilities, 
contained in Eqn. (Q), so does each column. Since P(Ai;Bj) > 0, we may set 
a 2 (Ai\ Bj) = P(Ai\ Bj); the a's are probability amplitudes. We now have 



[P 



AB 



a 



ABl 



( a 2 (A 1 ;B 1 ) a 2 (A i; B 2 ) 
a 2 (A 2 ;B 1 ) a 2 (A 2 ;B 2 ) 



a 2 (A i; B N ) \ 
a 2 {A 2 ;B 



(5) 



\« 2 (^ B i) ^ 2 (A N ;B 2 ) ... a 2 (A N ;B N ) J 
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As each row and each column sums to unity, this suggests that the a's 
are direction cosines, where A\, A2,...,An are N mutually perpendicular axes 
in TV-dimensional space and B\, B^t-.Bn is another set of orthogonal axes. 
This means that if another observable C with TV eigenvalues Ci, C?,,...Cm is 
introduced, the three sets of a's are connected through the relation 

N 

a(A i; C k ) =J2a(A i ;B j )a{B j ;C k ) (6) 
3=1 

which holds for direction cosines. Though we have a = VP, the most general 
case involves a phase, so that 

a(A t ; C k ) = i>{Ai-C k ) = ^P{A t ;C k )e^ (7) 

where ip is the phase. This choice also leaves the rows and columns of the matrix 
in Eqn. (^|) summing to unity. The two-way symmetry expressed in Eqn. ([!]) is 
now contained in the Hermiticity condition 

ip{A i ;C k ) = i,*{C k ;A i ). (8) 

Eqn. (||) now connects the probability amplitudes ip{Ai\ C k ), x(Ai> Bj) an d 
<f>(Bj\ C k ), and takes the form 

JV 

ip(Ai;C k ) =J2x(A i ;B j )ct>(B j ;C k ). (9) 

3=1 

We have set x{Ai\Bj) = xjj(Ai\Bj) and cj)(Bj;C k ) = ijj(Bj;Ck) to avoid 
confusion. This expression is evidently the well-known law of probability in- 
terference. It is the basis of the expansion of the wave function in terms of a 
complete set. To see this, we omit Ai, the label for the initial state and let 
Ck = x because the final eigenvalue spectrum is continuous. Then we set 

b j = X (A i ;B j ), o, oil 1:0,}. (10) 
We then obtain the familiar expression 

N 

iP(x) = Y / b j cf> j (x). (11) 

3=1 

2.2 Extension of Basic Results 

The following formulas are not explicitly given by Lande, but they are important 
to what subsequently follows. 

When A = C, Eqn. (0) becomes 

N N 

i>(A i ;A k )=Y,x(A i ;B j )<t>(B r ,A k ) = J2x(A i \B j )x*(B r ,A k ) = 5 ik (12) 

3=1 3 = 1 
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In the standard expression Eqn. flllD, the expansion coefficients bj are of course, 
probability amplitudes. However, the extra detail in Eqn. (^) shows that, far 
from being merely constants, these probability amplitudes x{Ai]Bj) have a 
structure. In fact, they are completely equivalent in character to the members 
of the basis set {<j)(Bj\ C k )}. Because of the fact that they are probability 
amplitudes, the expansion coefficients satisfy 

N N 

2lx(^;Bi)l 2 = EN a = 1 - ( 13 ) 

Usually, this result and the interpretation of the expansion coefficient as 
probability amplitudes are justified in a more roundabout way, while the Lande 
formalism justifies these results naturally and automatically. 

Since all the functions in the expansion Eqn. (^) are on the same footing, 
there are equivalent expressions for both xi^-i'i Bj) and 4>{Bj\ C k ). If we multiply 
the basic expansion Eqn. (^|) by x*{A-i\ Bi) and sum over i, we obtain 

N N 

J2x*(AilB l )^(A i ;C k ) = J2<t>(Br,C k )Y^X*(A i -,B l ) X (A i -,B j ). (14) 

i=l j = l i 

But in view of Eqn. (||) and Eqn. ( p^ ) , we obtain 

N 

<t>{B t \ C k ) = x(Bi;Ai)i>(Ai; C k ). (15) 



The other relation is found by multiplying both sides of Eqn. (||) by <p*(B m ; C k ) and 
summing over k. Due to the orthogonality of the functions <fi{Bj] C k ),we obtain 

JV 

x {A i; B m ) = i>( A i\ C k )<f>{C k ;B m ) (16) 
an expression that is better known in the form 

b m = J 4>* m (x)il)(x)dx. (17) 

3 THE EIGENVALUE EQUATION IN MATRIX 
FORM 

3.1 The Differential Eigenvalue Equation and Its Inter- 
pretation In The Lande Formalism 

We now proceed to an investigation of the connection between wave and ma- 
trix mechanics. We seek to obtain the matrix eigenvalue equation from the 
differential eigenvalue equation. 
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Though Lande does not explicitly treat this point, his formalism requires 
the continuous variable in the differential eigenvalue equation to represent the 
eigenvalues corresponding to the states being obtained upon measurement. The 
eigenvalue that appears as the label of the eigenfunction in the differential eigen- 
value equation describes the state that obtains before measurement. Thus, for 
example, in the single-particle time-independent Schroedinger equation 

ff(r)V fe (r) = E k ^ k {v) (18) 

E k represents the eigenvalue corresponding to the state that obtains before mea- 
surement, while r represents the result of the measurement. The eigenfunction 
ip k (r) is the probability amplitude that if the system is in a state corresponding 
to Ek, the value r is obtained upon measurement of the position. Similarly, in 
the equation 

d 

~ lh Q^^P W = P^PO ( X ) ( 19 ) 

Po represents the eigenvalue for the initial state, while x is the final-state eigen- 
value. Thus, in this equation, tp(x) gives the probability amplitude for measure- 
ment of the position of the system if the system is initially in a state character- 
ized by the momentum pq. 

In the Lande formalism, the state that obtains before measurement is well- 
defined and has been brought about by prior measurement. In a differential 
eigenvalue equation, the operator acts on the variables representing the final 
eigenvalue of measurement. Therefore, the differential eigenvalue equation is 
possible only if the eigenfunction has a continuous eigenvalue spectrum, so that 
the differential operator can act on a continuous variable. 

3.2 Prom The Differential to the Matrix Eigenvalue Equa- 
tion 

As is well known, the differential eigenvalue equation is transformed to a matrix 
eigenvalue equation in the following way. Let the eigenvalue equation be 

A(x)i)(x) = \i/)(x) (20) 

where A(x) is the operator and A is the eigenvalue. We expand ifj(x) in terms of 
the complete orthonormal set {(£>j(x)}, which has N non-degenerate members: 

N 

^(x)=J2Wx) (21) 
j=i 

where bj is a constant. We substitute into Eqn. (|2^); pre-multiplying by (j)* m , in- 
tegrating over x and using the orthogonality of the members of the set {4>j \ gives 

N 

Y,b j (cf> rn \A\cl> j )-\b m = Q, m = 1,2,3, ,N (22) 

i=i 
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With A m j = (0 m we obtain the following system of equations if we 

write out these equations for each value of m: 



[A 11 -X]b 1 -A 12 b 2 + ... +A 1N b N = 
A 2 ih + [A 22 - \}b 2 + ... + A 2N b N = 



(23) 



A N1 b 1 +A N2 b 2 + ... + [A NN - \}b N 
Hence the matrix eigenvalue equation is 



/ An — A A 12 
A 2 \ A 22 — A 



Nl 



liV2 



\ 



A 



NN 



( 6l \ 



\b N j 



= 



The vector 



[b] 



( 61 \ 

b 2 



\b N J 

is the matrix representation of the eigenfunction in the basis {<j>j}. 
The expectation value of the quantity R(x) is 

N N 

W = EE 6 * i2 *A- = M t [^] [6] 
t=l J=l 



where 



[i?] 



/ fill fil2 
fi21 fi22 



filJV \ 
R 2 N 



and 



\ Rni Rn 2 ••• fi^TV / 



R iJ = (<f> i \R(x)\<t> j ). 



(24) 



(25) 



(26) 



(27) 



(28) 



3.3 The Matrix Eigenvalue Equation in the Lande For- 
malism 

We now recast this procedure in terms of the Lande formalism. We start again 
with the differential eigenvalue equation: we write this out with as much infor- 
mation as possible: 
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A(x)ip(X k ; x) = \ k ip{\ k \x). 



(29) 



Here x is a general quantity which represents the final eigenvalue, whose spec- 
trum is continuous. We use Eqn. (pi) to express ijj(^k, x) as an expansion: 



N 



ip{\ k ;x) = Vx(Afc;S J )0(S 3 ;a;). 



3=1 



(30) 



We substitute this in the eigenvalue Eqn. fl29|). Then we pre- multiply by 
4>*(B m ,x) and integrate over x to obtain 



^ / <j)*{B m ;x)A(x)(j)(B j ;x)dxx{\k;B j ) = ^ A fe / (f>* (B m ; x)4>(B j ; x)dxx(\ k ; Bj). 
3=1 J 3=1 ^ 

(31) 

Due to the orthogonality of the probability amplitudes tf>, we obtain 



N 



^A mjX {\k;B j )-\ kX (\k;B j ) = 0, form =1,2, 3, AT (32) 

3=1 



where 



A mj = (<p(B m ;x) \A(x)\ (t>{Bj)x)) . 



(33) 



By writing down the equations for each value of m, we are led to the matrix 
eigenvalue equation 



/ An -Afe A 12 
A21 A22 — Afe 



\ Ajsix An2 
If we denote the eigenvector by 



Ain \ ( x(A fc ,Si) \ 
A2N 



Ann — Afe / 



x(Afc,S 2 ) 
V x(^k,B N ) J 



[C(Afe)] 



/ x(Afe,Si) \ 
x(Afe,S 2 ) 



V X (X k ,B N ) J 
then we obtain the matrix eigenvalue equation 

[A]K(Afe)]=Afe[£(Afe)]. 
The expectation value of the quantity R(x) is given by 



(34) 



(35) 



(36) 
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N N 



i=i j=i 



where 



«ij = <0(-B,;a;)|i?(x)|0(B,;x)) 



(38) 



/ x(A fc ;,Bi) \ 
x(A fe ;B 2 ) 



V x(A fc ;Sjv) / 



(39) 



and 



[R] 



( Rn R12 
R21 R22 



RlN \ 

R2N 



(40) 



y Rjsri Rn2 ■■■ Rnn J 

We remark that, for convenience, we have made a notational distinction between 
the arbitrary state [tp] and the eigenvectors of the operator [A] , which we denote 

JT-i A-i 

The eigenvalue equation (|3J) is the counterpart of Eqn. (|24|) but it contains 
extra detail which permits us to give a better interpretation of the quantities 
that appear in it. The elements of the eigenvector are the probability amplitudes 
x(Afc; Bj). In general, these probability amplitudes describe measurements from 
the initial state corresponding to the eigenvalue Afc to intermediate states which 
are different from the ones appearing as the final states in the differential eigen- 
value equation. Thus, the matrix mechanics description involves a complete 
set of intermediate states corresponding to an intermediate observable B. This 
observable may be chosen in a variety of ways; if it is chosen to equal the initial 
observable A, we obtain the results corresponding to using the eigenrepresenta- 
tion. 

The fact that the elements of the eigenvector are probability amplitudes 
is of course standard knowledge; however, we see here that these probability 
amplitudes have a structure and are far from being mere constants. As a result 
of this structure, we are able to obtain the eigenvectors of complex matrices by 
inspection, as will be seen in later sections. 

The present approach affords a unified treatment of measurements on quan- 
tum systems. Whatever the measurement we wish to make, we begin by talking 
in terms of the probability amplitudes for the various possible outcomes of 
the measurement. These can then be used to make the transition to matrix- 
mechanics. Since all results of measurement can be characterized by the prob- 
abilities that are associated with obtaining the various values than can result 
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from the measurement, this approach should always be possible. Therefore it is 
very general. 

In this paper, this approach is used to treat spin. It permits us to derive 
the spin vectors and operators from a more fundamental level of analysis even 
though the spin cannot be expressed in terms of the position coordinates. 

4 SPECIALIZATION TO SPIN-1/2 SYSTEMS 

4.1 Spin-1/2 Systems 

Consider a spin-1/2 system whose spin er is known to be up with respect to the 
unit vector a. A measurement of the spin along a new unit vector c is made. 
At the same time, a measurement of the quantity R(cr ■ c) which is a function 
of the spin projection cr-c is made. In order to describe the measurements, 
we introduce the probability amplitudes for the measurements. These are the 
equivalent of the wave function, except that they are discrete. We shall assume 
that all spins are measured in units of h/2. 

We want to use the Lande expression Eqn. @ to expand each of the am- 
plitudes. The initial and final states are spin states. The intermediate state 
is also a spin state. This arises from the fact that spin measurements along 
different quantization directions have the status of different observables. The 
advantage of this circumstance is that though the initial, intermediate states 
and final states formally have differing statuses, they in practice are connected 
by probability amplitudes of the same form. 

Let the states be labelled by the magnetic quantum numbers. Since each such 
number is with respect to a quantization direction, a superscript indicating this 
direction is also needed. Thus, (+^ )( c )) is the probability amplitude 

for the spin of the system being found upon measurement to be up with respect 
to c, after the spin had previously been ascertained to be up with respect to 
a. There are three other probability amplitudes; these are ip({+^)^; (— |)^ c '), 

and with obvio us interpretations. We 

note that orthogonality of the spin states corresponding to one quantization 
direction means, for example, that 

^((+\) a) ; (-\f y ) = (+l) ay ) = o (4i) 

According to Eqn. (^), we have the following expansions: 

4>a+\)®;i+\)®) = x((+^) ( " ) ;(+^) ( ^ ) M+^) (S) ;(+^ ) ) 

+x((+l) (7) ; (-^M^)®; (+5)®) (42) 
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^(-fhi+l^) - x((-i) ft ;(+^) (S) )0((+^) (S) ;(+i)^) 

+x((-^ } ; (-^) (S) )0((-i) (£) ; (+1)®) (44) 

and 

+x((-i) (2) ;(-i) (e) )0((-i) (S) ;(-i) ( ^). (45) 

Here b is another unit vector along which we can measure the spin projec- 
tion. The basis functions for the first and third expansions are the probability 
amplitudes (+|)®) and (+±)©). For the second and 

fourth expansions, they are </>((+|y b) ; (~|)®) and <p{{-\) {h) ; (-|) ®). 

According to Eqns. ( |30| ) and (p9|), the matrix representations of these prob- 
ability amplitudes -0 are 

= m(+b^(--^ ] = ( x((+ ^ (a);(+ ^ )(b)) 

and 



T ly 2' " 2' ' v 2' ,J \ v(7+iV a )- f-i)( b )) 



2 ' 



(46) 



M(-2) w ;(+2) w )] = [V'((-§) w ;(-2) w )] 



x((-i) w i(-j) (bl ) 



(47) 

We observe that provided the initial state is the same, the matrix state is the 
same. 

We now denote the quantum numbers as mi = +\ and TO2 = —\. The 

general probability amplitude ip may therefore be denoted by 0(m^;TOn ), 
where (i,n = 1,2). The same kind of notation will apply to the probability 
amplitudes x and <fi. 

Let the quantity R(cr ■ c) have the value r% when the spin projection is up 
(quantum number mi) with respect to the vector c, and T2 when it is down 
(quantum number 7712). Thus, the possible values of R are r n (n — 1,2). Sup- 
pose that the initial state corresponds to the quantum number rtii with respect 

to a. As the probability amplitude for obtaining r n is t/>(m^ , miif > ) 1 the expec- 
tation value of R is 

2 ~ ^2 

(R) = Y,m<\w$) r n . (48) 
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Since the expansions for -ip*(m[ a ^; m^) and ip(m\^; m^) are 



and 



it follows that 



where 



and 



= [^(m®;m®)]t[ii]^(mp;m®)] 



[R} = 



Ru R\2 



R21 R22 J ' 
with the elements of [R] being given by 

2 ^ 



Rjjl = Y,<P*{nif ) -m^)<p{ m p- m ^)r n 



n=l 

Written out explicitly, the elements of [R] are 

2 



Ru 



ri + 



'2 



+^((+^ (e) ;(-^®W(-^) (S) ;(-^)®)r- 2) 



«2i = ^(("^^K+^M+^^^+^Jri 

+^((-^) (e) ;(-i) ( ^)0((+i) (S) ;(-i) (?) )r2 



and 
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R 2 



<M(- 



4>((-^;(-^) 



r 2 . 



(58) 



When R = cr ■ c, then it is the component of the spin in the direction c. 
Hence, [R] is the matrix form of the operator for the component of the spin 



along the axis defined by c. 
square of the total spin. 



When R = cr 2 , then [R] is the matrix form of the 



4.2 Summary of the Various Possible Choices for the Ref- 
erence Vectors 

The results given above are the most general possible. We now consider several 
different cases that arise when we choose the vectors b and c in specific ways. 
Case (a): b ^ a and c ^ a. 

This is the most general case. In this case, the matrix representations are 



r " r 2' Vx((+I) (a) ;(-I) (b) ) 



and 

:(+-)&)} = [</>((--)£)• (--)% = i *(H) (a) ; (+3) (b) ) 
m\ 2 > ,k+ 2 ) n my 2 > ,y 2 > n y x ^_i^.^_i^ 

while the operator [R] has its elements given by Eqns. (pq ) - (§)■ 
Case (b): b = a. 

For this case, the matrix representations are 



(59) 



(60) 



M(+\)®;(+l)®)] = [iKC+^j = ( J ) (6i) 

(+lf y )} = (-^)®)] = ( J ) (62) 



and 



while the elements of [R] are Eqns. (|55|) - (|5§|) with b replaced by a. 
Case (c): b = c. 

The matrix representations are 

W( + i)® ;(+ i)®)] = W( + i)« ; (-i)«)] = ( $tft£j)ffl) ) 

(63) 

and 
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while the operator is 



Case (d): c = a. 

The matrix representations are 



(64) 



^ ( r 2 )' 



2 2 1 ' 2' ' 2' ! \ X ((-|)(«);(-& b )) 



(66) 



and 



w(4)«;(4)«)]=w(4)« ; (-i)«)]=( *«:p;<:pj ) 

(67) 

while the elements of [R] are given by Eqns. ( |55| ) - ( |58| ) with c replaced by a. 
We observe that despite the generalized form of the matrix representations, the 
probability amplitudes themselves are 

^>((+^; (+^f y ) = (-^f ] = i (68) 

and 

VX(+^) ft ; = iW-^)®; (+^) ( * } ) = o. (69) 

Case (e): b = a and c = a. 
The matrix representations are 

[V>((+^ A (+-)% = MC+^s (-i)^)] = ( J ) (70) 
(+^) ( ^)] = (-^)% = ( J ) (7i) 



and 



while the operator is 
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We observe that Cases (a)-(c) are just different ways of treating exactly 
the same measurements and computing the same expectation value. The initial 
projection axis and the same final projection axis are all the same; the situation 
corresponds to the general case. On the other hand, Cases (d) and (e), which 
are also equivalent, both refer to a special case. Since a = c, both describe a 
repetition of the spin projection measurement along a given axis. 



4.3 Derivation of the Pauli Spin Matrices and Vectors 

The expressions for [R] allow us to deduce the matrix forms of various operators. 
Thus if R is the spin projection itself, then r\ = 1 and r 2 = —1, so that from 
Case (e), we are able to deduce the operator for the spin projection along the 
axis of quantization as 



1 

-1 



(73) 



Eqn. @ is of course the Pauli matrix for the z component of the spin. 

Even though the z axis has not been explicitly introduced, we are using z 
as a subscript for a good reason. The generalized form of this operator, to be 
derived below, will be defined with respect to an arbitrary unit vector, say c. It 
will be seen that in an appropriate limit, this generalized operator will reduce to 
the Pauli matrix Eqn. (k3h. This circumstance justifies our use of the subscript 



Now, since the eigenvalue of a 2 is 3, we have r\ = r-2 = 3. Therefore the 
operator for the square of the spin is given by Case (e) as 

M -(11) (74 > 

We emphasize that these operators apply when b = c and when a = b = c. 

Once we have the z component and its eigenvectors, we invoke standard 
theory to obtain the x and y components. Thus, we derive the matrix forms 
of the operators for the x and y components of spin by first obtaining the 
ladder operators [er+] and [cr_] . These are deduced from their actions on the 

eigenvectors of [crj. The eigenvectors of [a z ] are [£ + ] ~ ^ q ^ corresponding 

to the eigenvalue +1, and [£_] = ( ^ ^corresponding to the eigenvalue —1. 

In terms of the ladder operators, the operators for the x and y components 
of the spin are 

[<r*] = \(W+] + [*-]) (75) 

and 

Wv\ = -^W+] -[*-])■ (76) 
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The effects of [er+] on the eigenvectors of [<r 2 ] are 

["+][£+]= (77) 

and 

[* + }{u = m+i (78) 

Using these properties, we find that 

M=2(° I). (79) 

Similarly, using 

[<7_] IU = (80) 

and 

[<?-] IU = 2[f + ] (81) 

we find that 

M=2(5 °). (82) 

Hence 

(83) 
and 

w=(; -;). (84) 

These are the Pauli spin matrices, which we have here derived from first 
principles. We see that the Pauli spin vectors and matrices arise as a special 
case of spin measurements. They result if we make a second measurement of 
spin without changing the reference direction, and while using the initial/final 
direction to define the intermediate states which appear in the expansion of the 
relevant probability amplitudes by Eqn. (^). 

We need to obtain the most general possible vectors and operators for the 
description of spin. In order to obtain these quantities, we need to obtain explicit 
expressions for the probability amplitudes \ an d 4>- 




4.4 Explicit Formulas for the Probability Amplitudes 

In the last section, we have obtained the matrix operator for the spin, for the 
special case corresponding to the Pauli vectors and matrices; this operator is 



1 \ -/ -i \ t 1 

1 +J h o +k -1 



(85) 
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Consider the operator obtained from a by taking its dot product with a 
unit vector a. A measurement of spin along the new direction defined by the 
unit vector a gives the values ±1. Let the polar angles of a be {9,(p), so that 
a = (sin 9 cos ip, sin sin ip, cos 9) . The component of the spin along a is <x • a. In 
view of Eqn. ( p^ ) and of the components of a, the matrix form of this operator 
is 

r ~ / cos 9 smOe~ %v \ , . 

^■^ = {sm9e^ -cos 9 J' (86) 

The eigenvalues of this operator are +1 with eigenvector 

os (9/2 
'sin 9/2 



= 1 j: s JJL i m 



and —1 with eigenvector 



-e lv cos6»/2 

From the matrix eigenvalue equation (|34|), we deduce that the elements of 
the eigenvectors are of the form x(^k',Bj), where Afc is an eigenvalue of the 
matrix, and corresponds to the state that precedes measurement, while Bj is an 
eigenvalue corresponding to another state. Thus, these probability amplitudes 
have the structure of the x(^k] Bj). 

In general B is not the observable being measured and is also different from 
the observable corresponding to the initial state. But for the case of spin mea- 
surements, provided that B represents spin projection measurements along a 
new direction, it has the status of another observable. We therefore conclude 



that in the vectors [^-J and [£ 2 ]> Eqns. ( p7| ) and (88) respectively, the elements 
are probability amplitudes which refer to measurements from an initial state 
in which the spin is up or down with respect to a . In [^J, the elements are 
probability amplitudes corresponding to an initial state characterized by the 
spin projection being up with respect to a; the final states correspond to the 
eigenvalues +1 and —1 with respect to the unit vector b. In [£ 2 ], the elements 
are probability amplitudes corresponding to the initial state characterized by 
the spin projection being down with respect to a; the final states correspond to 
the eigenvalues +1 and —1 with respect to the unit vector b. Hence we deduce 
from Eqn. @ that 

x{{+\f ) -,{+\f ) )=™s9/2 (89) 

and 

X {{+\f ) :,{-\f ) )=^sin9/2. (90) 
Similarly, we deduce from Eqn. that 



x{{-\) {2) ;{+\f ] ) =sinfl/2. (91) 
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and 

x({~)^\(~)^) = ^0^6/2. (92) 

For the time being, we do not need to know what b is, though we shall 
deduce it later. As b may be a special vector, we cannot assume that the 
probability amplitudes Eqns. ( |S9| ) - ( p2[ ) are general. But armed with these 
special amplitudes, we are able to obtain general expressions for the probability 
amplitudes using Eqn. (|^). 

At this point, we note a simplification which arises from the fact that in 
all the probability amplitudes under consideration, the initial and final states 
both correspond to spin projections. This means that in reality, the probability 
amplitudes ip, x an d 4> have identical forms. 

Let c be a new unit vector defined by the polar angles (0' , ip'). Then if we 
make measurements of spin from c to the unknown quantization direction b, 
the probability amplitudes are 

x{{+\f ) ;{+\f ) )=^6'/2 (93) 

x{{+\f ) -A-\f ) )=e^ ^9'/2 (94) 

x((-^;(+~)^)= S m072 (95) 

and 

x((™)®;(-i)®) = -e*W/2. (96) 

We can now use Eqn. (|9|) to eliminate the unknown vector b, and thereby 
find the general probability amplitudes x- Using the Hermiticity condition Eqn. 
(p|), we find that 

x({+\f ) ;{+\f ) ) = x((+i) ft ;(+i) (S) M(+i) (S) ;(+^^) 

+x((+i) (2) ;(-i) (S) )x((-i) (S) ;(+i) ( ^) 



cos 



6>/2 cosO' /2 + e i{v - v,) sin 0/2 sin0'/2. (97) 



Similarly, we find that 



x((+^) ( ^ ) ;(-^) ( ° ) ) = cos0/2sin072-e l(¥, ~ v ' ) sin0/2cos072 (98) 



x ((__)(a).( + _)(c)) = sin0/2cos0y2-e l ^~^cos0/2sm0'/2 (99) 



and 
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x((-\)®;{-h®) = sin 9/2 sin 0'/2 + e^ - ^ cos 0/2 cos 6»'/2. (100) 



The expressions Eqns. ([97]) - (10C ) are the general expressions for the probability 
amplitudes which we seek. 



4.5 Verification of the General Expressions 

We need to check that the expressions for the probability amplitudes we have 
obtained are correct. If they are, then they should satisfy the Lande expansion 
Eqn. (j|). Let e be a new unit vector defined by the polar angles (6", ip"). Then 
according to Eqns. @ - fllO(i|) , 



x (( + i)(a) ; = cos 6»/2 cos 6"/2 + e^-^ sin 0/2 sin 6"/2 (101) 



x((+^) (2) ; = cos 6»/2 sin 6»"/2 - e^-?") sin6» /2cos6>" /2 (102) 



x ((_i)(S) ; = sin 6>/2 cos 6»" /2 - e^^"-" cos<9/2sin6>72 



3" /O — J{<P-<P") , 



(103) 



and 



XCI-^)^;!-^)^) =sin0 / 2sinr /2 + e i(v - v " ) cos6'/2cos6'"/2. (104) 



Also, 

X ((+I)(S>; (+!)£)) = cos6»"/2cos6 , 72 + ^'-v") sin 6»"/2 sin 6*72 (105) 



x((+^) (a) ; = cos^72sin&72- e*^""^ sin" 5/2 cos 0'/2 (106) 



x((-^)®;(+^) (S) ) = sin6>72cos6'72-e i ^"- v ' ) cos5"/2sin6'72 (107) 



x((-^) (S) ;(-^) ( " ) ) = sin e'72 sin 6*72 + e i(¥> "- v,) cos 6>'72 cos 6*72. (108) 
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Now, 



x«+\)®',(+\)®) = x((+i) ft ;(+i) (a) )x((+^) (a) ;(+i) ( ^) 

+x((+^) ft ; (-^)®)x((-i) (a) ; (+~)®) 

= cos 6»/2 cos 9' /2 + e i{ ^- <p ' ) sin 9/2 sin 9'/2. (109) 

W he n we substitute for the various probability amplitudes from the list Eqns. 
(101) - (108), we recover Eqn . (j97[). The same holds for the other probability 
amplitudes Eqns. ( p8|) - ( |l00| ). This confirms the consistency of the expressions 
for the probability amplitudes and the validity of the Lande expansion Eqn. (^). 



4.6 Expressions for the Probabilities 

The probabilities corresponding to the probability amplitudes Eqns. ( |97j ) - ( |100[ ) 

are 



P((+l)®; (+5)®) 



cos 2 (6>/2 - 9' / 2) - sine 1 sine 1 ' sin 2 (<p/2 - <p'/%).10) 



Pd+h^l = sin 2 (6»/2 - 0'/2) + sin(9sin0' sin 2 (<p/2 - y//2) (111) 



P((_i)(£) ; (+1)< C )) = sin 2 (6»/2 - 6»'/2) + sin(9sin0' sin 2 (<p/2 - <p'/2) (112) 



and 



P((_I)(S) ; (-l)W) = cos 2 (672 - 0'/2) - sin 9 sin(9' sin 2 (<p/2 - <//2). (113) 



2' ' v 2 

The cosine of the angle 8 between the vectors a and c is 

cos 8 = a ■ c = cos(6 i — 9') — 2 sin 9 sin 9' sin 2 (</? — <p')/2. 

Hence, 

cos 2 8/2 = cos 2 (6»/2 - 9'/2) - sin 9 sin 9' sin 2 (<p/2 - tp'/2) 
We see therefore that 

p((+\) & ; (+\f y )) = P((-\) a) ; (-\f } ) = cos 2 e/2. 

In the same fashion, we deduce that 

P((+\) ah , (~lf ] ) = P((-l) {2 \ (+5)®) = sin 2 8/2. 



(114) 



(115) 



(116) 



(117) 
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As we expect, 

P((±i)® + P((±\)®; (-i)®) - l. (118) 

Thus, the probabilities are normalized to unity as we expect. 

We are now in a position to deduce the vector b. It is straightforward to see 
that b is the vector with angles 0=0 and tp' — ir. We expect therefore that 
if we make these substitutions in the generalized probability amplitudes Eqns. 
( p7| ) - ( |100| ), we should recover the specialized amplitudes Eqns. ( ]8^ ) - (B^). 
This is precisely the case. 



4.7 Explicit Formulas For the Operators 

We can now obtain explicit general formulas for the various operators. We want 
to obtain the explicit expressions Eqns. ( p5| ) - (|5^) which are the elements of the 
general operator [R] . Since these are given in terms of vectors b and c, we have 
to redefine the angles so that they refer to these vectors. We define the unit 
vector b by the angles (0,tp), while the unit vector c is defined by the angles 
(9' , tp'). It is important to realize that the vector b in this subsection is different 
from the vector b in Subsections 4.4 and 4.6. The vector in Subsections 4.4 and 
4.6 is a particular instance of the general vector in this section. 

Using the fact that the (fi's and x' s have identical forms, we find that 



<M(+J) (S) ;(+7j) ( ^) = cos6»/2cos072 + e' i(¥ '- ¥ '' ) sin6l/2sin072, (119) 



(__)©) _ cos6»/2sin6'72-e 4(¥, - ¥, ')sin6'/2cos072, (120) 



= sin 6V2 cos 6*72 -e 1 ^-^ cos (9/2 sin 6>72 (121) 



and 



- sin 6>/2 sin 6»72 + e !(v - v,) cos 0/2 cos 6*72. (122) 



Hence, using Eqns. ( J55[ ) 
elements are 



we find that for the general operator [R] the 



Rn 



<H(- 

[cos 2 



9')/2 - sin 9 sin 6' sin 2 (ip - cp') /2}n 
-9')/2 + sin 8 sin 6' sin 2 (tp - tp') /2]r 2 



(123) 
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R\2= (y ^ 1 ^ T2 ^j (sin 9 cos 9 1 — sin 9' [cos 9 cos(p — <p') + isin(Lp — <p')]) (124) 

i?2i = ^~ ~2 — (sin ^ cos 6 1 ' — sin ^' [cos 6 1 cos(^ — (jC 7 ) — z sin(cp — 9^')] ) (125) 
and 

i?22 = [sin 2 (6»-6»')/2 + sin6'sin6l'sin 2 (^-^ / )/2]ri 

+ [cos 2 (0 - 0')/2 - sin6» sin 9' sin 2 (<p - tp')/2]r 2 - (126) 

4.8 Generalized Form of [a z ] and of its Eigenvectors 

Certain special operators are now readily obtained. If we seek the operator for 
the projection of the spin along the new axis defined by c, then we set r% = 1 
and r2 = —1. Hence the operator for the component of spin in the direction c is 



(<r~)n (cr-)i2 
(&~hi (0^)22 



(127) 



where 



(cr-) u = cos(6> - 6^-2 sin 6> sin 6>' sin 2 (<^-</)/2 
= cos 9 cos 9 + sin 9 sin 9' cos (<p — <p') 

(0^)12 = sin^cos^' — sin 0'[cos 9 cos(p — ip') + ism(p — p')\ 
(0^)21 = sin 9 cos 6 — sin 6*' [cos 9 cos (95 — tp') — ism(p — tp')] 



and 



(0^)22 



9' — sin0sin0'cos(<£ — ip'). 



(128) 

(129) 
(130) 

(131) 



The operator for the square of the spin is obtained from the general formula 
by setting n = r 2 = 3. Hence 



\* 2 \ = 



3 
3 



(132) 



The operator defined by the elements Eqns. (128) - (131) is the most general 
for the component of the angular momentum along the axis of quantization. We 
need the operators for the x and y components. We can obtain these through 
the raising and lowering operators. But in order to deduce these operators, we 



need the eigenvectors of the operator Eqn. (127), whose elements are given by 
Eqns. ( |128[ ) - (131). We could obtain these by direct calculation, but, given 
the form of the operator, this is likely to be tedious. However, we can deduce 
these eigenvectors by arguments based on the interpretation of the eigenvalue 
equation (p3). 
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The elements of the eigenvectors are probability amplitudes. From the ma- 
trix eigenvalue equation (Q), we see that their form is given by Eqns. (|97]) - 
( p-OOp . We only need to identify the initial and final directions corresponding 
to these probability amplitudes. This we can do by inspecting the five cases 
contained in Eqns. ( |56| ) - (f72|). 

The generalized operator [cr^], Eqns. 

two unit vectors b and c whose polar angles are (9, ip) and (9 1 , ip') respectively. 
Hence the elements of its eigenvectors must also be functions of these four angles 
only. The general case when the matrix and the eigenvectors are functions of 
two vectors only is possible if a = c, but b 7^ a. This corresponds to Case 
(d) among the possibilities Eqns. ( j58| ) - ([72]) ■ The expectation value of the 
quantity R is in this case equal to the eigenvalue of R corresponding to the 
initial state. This is true because c = a means that the measurement is a 
repeat measurement. But in order for the matrix expectation value to yield 
this result, it must be that the vector is an eigenvector of the operator. We 
therefore conclude that the eigenvectors contain as their elements probability 
amplitudes with initial states corresponding to the unit vector a and final states 
corresponding to the unit vector b. As a = c = (sin 8' cos <p', sin 9' sin ip', cos 9'), 
and b = (sin 9 cos ip, sin9 sin ip, cos 9), we obtain the required eigenvectors as 



128) - fll3ip, is a function of only the 



and 



x(H) (a) ;H) (b) ) 
x((+i) (2) ;H) (b) ) 



cos 9' jl cos 9/2 + e^''^ sin 6' / '2 sin (9/2 
cos 0'/2 sin 0/2 - e 



*(¥>'- 



sin 6»72 cos 0/2 
(133) 



, (-) / x((-!) (a) ;(+|) (b) ) \ _ ( sin 8' /2 cos 9/2-e l ^'~^ cos 6*72 sin 9/2 
x ((_i)(a). (-i)W) J~\ sin8' /2sin9/2 + e i ^'-^ cos cos 0/2 

(134) 

When we test these eigenvectors with their operator, we find that indeed the 
eigenvalue equations 

c c c 

and 

Wc]l^ ) ] = (- 1 )i^ ) ] (136) 

are satisfied. 



4.9 Generalized Forms of [a x ] and [a y ] 

We are now able to deduce the generalized forms of the ladder operators and 
hence of [a x ] and [cr y ]. For this purpose we use Eqns. (77), (J7§), ( p0| ) and ( pi] ) 
with the eigenvectors of [ct^], Eqns. ( |133] ) and (134). We find that the elements 
of [er + ] are 



( (J +)n = 



sm f cos 



?' cos(tp — <p') + cos ( 



sm( 



+ i sin 9 sin((p' — ip) (137) 
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(>+)l2 = 



2 [cos 2 0/2 cos 2 6'/2 + sin 2 0'/2 sin 2 0/2] cos(v5 - p') + sin sin 0' 
+2i[sin 2 6/2 sin 2 6>'/2 - cos 2 0'/2 cos 2 0/2] sin(<p' - (138) 



(c7+) 21 = -2[sin 2 0/2cos 2 072 + sin 2 072cos 2 0/2]cos(c / 3-</)+sin0sin0' 



/21 

-2i[sin 2 0'/2 cos 2 0/2 - sin 2 0/2 cos 2 0'/2] sin(p' - (139) 
and 

((T + ) 22 = sin0cos0'cos(<p — <^') — cos0sin0' — i sin sin(y/ — y>). (140) 
The elements of [<r_] are 

(<7_) n = — sin cos 0' cos(<^' — ip) + cos sin 0' — isin0 sin(<^' — ip) (141) 



(ct_) 12 = -2[cos 2 0/2sin 2 072 + cos 2 072sin 2 0/2]cos(^'-^)+sin0sin0' 
-2i[cos 2 072 sin 2 0/2 - cos 2 0/2 sin 2 072] am(<p' - ip) (142) 



(ct_) 21 = 2[sin 2 0/2sin 2 072 + cos 2 072cos 2 0/2]cos(^' - ip) +sin0sin0' 
+2?;[cos 2 0/2 cos 2 072 - sin 2 072 sin 2 0/2] sin(<p' - ip) (143) 

and 

(cr_) 22 = sin0cos0'cos((^' — ip) — cos sin 0' + isin0sin(<// — <p). (144) 

Therefore the generalized operators for the x and y components of the spin 
have the elements 

{ (J x)n = — sin cos 0' cos (tp' — ip) +sin0'cos0 (145) 
(cx)i2 = cos0cos0'cos(<// — ip) + sin sin 0' — i cos0' sm(p' — p) (146) 



(f x)2i — cos cos 0' cos(<// — <£>) + sin sin 9' + i cos 0' sin(<^' — p) (147) 

(o" x ) 22 = sin cos 0' cos{<p' — p) — sin 0' cos (148) 

{&y)ii = sin0sin((^' — ip) (149) 

(cj/)i2 = —icos(tp' — ip) — cos0sin(<^' — ip) (150) 
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(cy)2i = i cos(tp' — <p) — cos 9 sin(<^' — ip) (151) 

and 

(o"w)22 = -sin#sin((p' - y>). (152) 

We verify that 

[a x f = [o y Y = [^] 2 = 1. (153) 
Also, the matrices anticommute: 

[[a x ], [a y ]} + = [[a y ], [a~)} + = [[a-], [a x }]+ = 0. (154) 
The commutators are 

[M, [<t v ]] = 2i[a~] (155) 

[[<r v ], [a~]} = 2i[a x ] (156) 

and 

[[a-],[a x }}=2i[a y ]. (157) 

Thus the generalized spin matrices obey the relations which the well-known 
specialized forms, the Pauli spin matrices, Eqns. ([73]), ( |S3| ) and (Q) obey. 

At this point, it is worth observing that the x and y components of the 
spin have been so identified purely because they are connected to the ladder 
operators by the standard relations Eqns. ( f75| ) and (|76|). We expect that in the 
limit in which the generalized form [<r— ] reduces to the Pauli spin matrix [cr z ], 
these operators will also reduce to the corresponding Pauli forms. 



4.10 Special Forms of the Vectors and Operators 

We need to verify that in the appropriate limits, the generalized results obtained 
here reduce to the ones corresponding to the special cases. Indeed, we find that 
we obtain the Pauli spin matrices by setting 9 — 9' and ip = ip' . 

We obtain the results Eqns. ( ]86| ) - ( p8{ ) by keeping a untouched, but setting 
9 = and tp = n. We then find that 



sin 9'e lv ' - cos ( 



(158) 



sin6»'/2 



and 
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We see that we do indeed recover the operator Eqn. (|86|) and the eigenvectors 
Eqn. @ and Eqn. (||). So the results given by Eqns. ( |l58| ) - ( |l60| ) are 
consistent. 

In the same limit, the ladder operators change as follows: 



[*+] 



and 



Also 



while 



sin#' 
2 sin 2 e'/2e lv ' 

sin 8' 
-2 cos 2 6'/2e iv ' 



-2 cos 2 e'/2e 

- sin 9' 

2 sin 2 0'/2e- 

— sin 9' 

-p-i<p' ms ft' 



■lip 



■lip 



-lip 



-te 



lip 







(161) 



(162) 



(163) 



(164) 



In order to obtain the Pauli spin matrices and vectors from the general 
formulas we have had to set 8 = 8' ,<p — ip' . The operators in Eqns. ( p3[ ) and 
( [34} ) appear from the general cases in the same limit as the Pauli operators and 
vectors. Therefore in order to get them from the operators in Eqns. (163) and 
(164), we have to set 6' = 0, and ip' = ir, since we have set 9 = and ip = tt to 



get Eqns. (163) and (164). This is easily confirmed. 



4.11 Connection with Standard Results 

We have now obtained the most general forms of the vectors and operators for 
spin-1/2 systems. It is important to clarify the relation between these quantities 
and the standard ones. It is of particular interest to expand on how they are 
connected to the Pauli matrices and vectors. 

The Pauli spin matri ces E qn s. j73l ), (|83] ) and (|84|) are particular c ases of 
the general cases Eqns. (|l28| ) - (|l3l|), Eqns. jl45|) - ( fus) ) and Eqns. (|l49| ) - 
( |l52|) . As we have seen, they arise from the generalized results when 9 = 9' and 
<P = V' ■ 

Now, the generalized operators correspond to the probability amplitudes 
for spin measurements having been expanded by the Lande formula using an 
arbitrary vector to supply the intermediate states that appear in the expansion. 
The less generalized formulas Eqns. (^) - (^||) use the specialized vector d 
whose polar angles are 8 = and (p = tt to define the intermediate states. The 
Pauli quantites result when the arbitrary vector used in the expansion equals 
the direction of initial quantization. 

The probability amplitudes which constitute the elements of the vectors 
which comprise states in matrix mechanics are not just constants but have a 
structure. This means that in the Pauli spin vectors 



1 



and [£_] 
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the constants 1 and which appear as the elements of the vectors are probability 
amplitudes. When the probability amplitude is unity, it represents a repeat 
measurement; when it is zero, it represents a measurement to a state orthogonal 
to the initial state. Thus from [£,] and [£_], we see that 

X ((±i)£) ; (±I)(a) )= l 

while 

X ((±I)£);( T I)£))= . 

Here a is an arbitrary direction vector. 

It is clear from Cases (a) - (c) enumerated above that the problem of de- 
scribing spin measurements can be approached in a variety of ways. The task of 
computing the expectation value of a quantity depending on the spin projection 
can be achieved in several ways. The most general way is illustrated in Case (a) , 
where the most general vectors and operators are used. But the same results 
are obtained if we let the intermediate direction equal the initial. Then the 
vectors simplify to the well-known "up" and "down" vectors, and the fact that 
the measurement is from one direction to another is contained in the operator, 
which has its most general form. Another way of simplifying the calculation 
is to let the intermediate direction equal the final direction. In that case, the 
operator simplifies to the Pauli matrix for the z component of spin, while the 
vectors remain general. 

In the special case of computing the expectation value when the measurement 
refers to the same initial and final direction, there are two possibilities. The 
initial and the final directions are the same. If the intermediate direction is 
not equal to the initial/final direction, then the vectors remain general and the 
operator also remains general. This is the hard way of doing it, because if the 
intermediate direction is chosen to equal the initial/final direction, then the 
vectors reduce to the "up" and "down" vectors, and the operator to the Pauli 
matrix for the z component of spin. 

We observe that when we choose the intermediate direction to equal the 
initial direction, as in Cases (b) and (e), this is formally equivalent to expanding 
the probability amplitudes in the eigenrepresentation, because this is when the 
operator matrices are diagonal. 

5 SUMMARY AND DISCUSSION 

In this paper we have used the Lande interpretation of quantum mechanics to 
investigate the connection between wave and matrix mechanics. We have used 
Eqn. (|^) to obtain a more detailed expression for the expansion coefficients that 
appear in the usual expansion of a wave function in terms of a basis set. We 
have reversed the expansion to obtain an expansion of the basis functions in 
terms of the expansion coefficients and the wave functions. 
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We have studied the transition from wave to matrix mechanics within the 
framework of the Lande interpretation of quantum mechanics. By writing down 
the matrix eigenvalue equation in a fully detailed way, we have demonstrated 
that the elements of the eigenvectors resulting from the solution of the eigen- 
value equation are probability amplitudes with a structure. These probability 
amplitudes describe a measurement starting from the state characteristic of the 
eigenvalue of the matrix to one of the states corresponding to the intermedi- 
ate observable brought into operation by expanding probability amplitudes by 
means of Eqn. @. This structure has been exploited to obtain the eigenvectors 
of complex matrices by inspection. 

In the derivation of the matrix treatment of orbital angular momentum, we 
normally use the spherical harmonics as basis functions in order to obtain the 
matrix angular momentum states and the matrix elements of the operators. On 
the other hand, spin does not have a treatment proceeding from the solution 
of a differential equation because. But we have shown that by using the prob- 
ability amplitudes for spin projection measurements, we can obtain the matrix 
treatment of spin in the same way that we obtain the matrix forms of the or- 
bital angular momentum quantities. Thus, we have derived from first principles 
the vectors and operators for spin 1/2. With only the assumption that the 
matrix operators and vectors thereby resulting have the properties of angular 
momentum operators and states, we have derived explicit expressions for the 
probability amplitudes. Moreover, we have derived the most general forms for 
the spin operators. We have not seen these forms in the literature; at any rate, 
this method of deriving these quantities is new. In the appropriate limit, our 
operators and vectors reduce to the Pauli spin vectors and matrices, which are 
thereby seen to be special forms of the generalized results we have presented. 

We make an interesting observation. A matrix-mechanics vector is charac- 
teristic only of the eigenvalue corresponding to the state that obtains before 
measurement, irrespective of the final state. Thus even zero probability am- 
plitudes may be represented by non-zero vectors. An example is afforded by 
Eqns. ( p8| ) and (^). The matrix state is a compact way of expressing discrete 
probability amplitudes, since a particular vector represents all those probability 
amplitudes which belong to the same initial state. 

Apart from the new forms of the matrix quantities, we have devised a general 
method for obtaining the generalized forms of spin operators and vectors. It 
is natural to seek to extend this approach to the derivation of vectors and 
matrix operators for other values of spin. We have successfully applied the ideas 
presented here to the case J = 1. In this fashion, we have obtained for J = 1 
quantities more generalized quantities than any so far seen in the literature. 
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